Under consideration for publication in J. Fluid Mech. 



1 



Anisotropy and cyclone-anticyclone 
asymmetry in decaying rotating turbulence 

By F. MOISY 1 , C. MORIZE 1 , M. RABAUD 1 , 
J. SOMMERI A 2 

1 Universite Paris-Sud 11, Universite Pierre et Marie Curie, CNRS; Laboratoire FAST, Bat. 

502, F-91405 Orsay, France 
2 Coriolis/LEGI, 21 avenue des Martyrs, F-38 000 Grenoble, France 

(Received 14 September 2009) 

The effect of a background rotation on the decay of homogeneous turbulence produced 
by a grid is experimentally investigated. Experiments have been performed in a channel 
mounted in the large-scale 'Coriolis' rotating platform, and measurements have been 
carried out in the planes normal and parallel to the rotation axis using particle image 
velocimetry. After a short period of about 0.4 tank rotation where the energy decays 
as £ -6 / 5 , as in classical isotropic turbulence, the energy follows a shallower decay law 
compatible with i -3 / 5 , as dimensionally expected for energy transfers governed by the 
linear timescale f2 _1 . The crossover occurs at a Rossby number Ro ~ 0.25, without 
noticeable dependence with the grid Rossby number. After this transition, anisotropy 
develops in the form of vertical layers where the initial vertical velocity remains trapped. 
These layers of nearly constant vertical velocity become thinner as they are advected 
and stretched by the large-scale horizontal flow, producing significant horizontal gradient 
of vertical velocity which eventually become unstable. After the Ro ~ 0.25 transition, 
the vertical vorticity field first develops a cyclone-anticyclone asymmetry, reproducing 
the growth law of the vorticity skewness, S u (t) ~ (fit) , reported by Morize, Moisy 
& Rabaud [Phys. Fluids 17 (9), 095105 (2005)]. At larger time, however, the vorticity 
skewness decreases and eventually returns to zero. The present results indicate that the 
shear instability of the vertical layers contribute significantly to the re-symmetrisation of 
the vertical vorticity at large time, by re-injecting vorticity fluctuations of random sign 
at small scales. These results emphasize the importance of the initial conditions in the 
decay of rotating turbulence. 



1. Introduction 

Turbulence subjected to solid body rotation is a problem of first importance for engi- 
neering, geophysical and astrophysical flows. Its dynamics is dictated by a competition 
between linear and non-linear effects. Linear effects, driven by the Coriolis force, include 
anisotropic propagation of energy by inertial waves, preferentially along the rotation axis 
(hereafter called 'vertical' axis by convention), on the timescale of the system rotation 
fi -1 (Greenspan 1968). Nonlinear interactions, on the other hand, are responsible for en- 
ergy transfers towards 'horizontal' modes (Jacquin et al. 1990; Waleffe 1993). For infinite 
rotation rate, i.e. for vanishing Rossby number, these inertial waves reduce to Taylor- 
Proudman columns, corresponding to a two-dimensional flow invariant along the rotation 
axis. Importantly, this two-dimensional flow is not two-component in general, because 
the third (vertical) velocity component, insensitive to the Coriolis force, behaves as a 
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passive scalar field transported by the horizontal flow. This 'passive' vertical velocity, 
originating from the initial conditions in the case of decaying turbulence, may however 
become 'active' through shear instabilities at small scale. This mechanism may have con- 
siderable importance in the nature of the decay and the partial two-dimensionalisation 
of an initially 3D turbulence subjected to background rotation. 

The present paper reports an experimental study of the influence of the background 
rotation on the decay of an initially isotropic turbulence. Turbulence is generated by 
translating a grid in a channel mounted on the large-scale 'Coriolis' rotating platform. 
The aim of this paper is first to characterize in details the decay law of the energy, in 
a situation where the effects of the lateral confinement can be neglected. This situation 
contrasts with the previous experiments by Morize, Moisy & Rabaud (2005) and Morize 
& Moisy (2006a), performed in a rotating tank with an aspect ratio of order 1, showing 
significant confinement effects. Second, the anisotropy growth is investigated, with the 
aim to characterize the influence at large time of the initial vertical fluctuations on the 
vertical vorticity statistics. 

For turbulence subjected to moderate rotation (Rossby number Ro ~ 0(1)), the lin- 
ear and nonlinear time scales are of the same order, resulting in a complex interplay 
between linear energy propagation by inertial waves and anisotropic energy transfers by 
nonlinear interactions. This complexity is unavoidable in decaying rotating turbulence 
starting at large Rossby number, in which the instantaneous Rossby number decays and 
crosses O(l) at some transition time. At this time the effects of the rotation, namely the 
anisotropy growth and the cyclone-anticyclone symmetry breaking, become significant 
and accumulate as time proceeds. Accordingly, the statistical properties of rotating tur- 
bulence at large time are the result of the turbulence history integrated from the initial 
state, and may therefore depend on the details of the initial state. Generic properties 
should however be expected if the initial state is 3D isotropic turbulence with Ro 3> 1, 
which is the situation examined in the present paper. 

Because of the fast growth of the vertical correlation due to inertial wave propaga- 
tion (Jacquin et al. 1990; Squires et al. 1994), confinement along the vertical axis plays 
a significant role in the dynamics of rotating turbulence. This vertical confinement is 
usually present either through solid or free-slip boundaries in experiments, or through 
periodic boundary conditions in numerical simulations. In all cases, comparisons with 
homogeneous turbulence in idealized unbounded systems should be made carefully. One 
consequence of the vertical confinement by solid or free-slip boundaries is to preferentially 
align the axis of the vortices normal to the walls, therefore reinforcing the two-dimensional 
nature of the large scales, as observed by Hopfinger, Browand & Gagne (1982) and Gode- 
ferd & Lollini (1999). Second, confinement selects a set of discrete resonant inertial modes 
(Bewley et al. 2007), which may couple to the small-scale turbulence. Third, an extra 
mechanism of dissipation of the inertial waves take place in the boundary layers, acting 
on the Ekman time scale /i^fi) -1 / 2 , where h is the confinement scale along the rotation 
axis, which may dominate the energy decay at large time (Phillips 1963; Ibbetson & 
Tritton 1982; Morize & Moisy 2006a). 

The most remarkable feature of rotating turbulence is the spontaneous emergence of 
long-lived columnar vortices aligned with the rotation axis (Hopfinger et al. 1982; Smith 
& Waleffe 1999; Longhctto et al. 2002). Both linear and non- linear mechanisms have 
been proposed to explain the formation of those columnar structures (Cambon & Scott 
1999; Cambon 2001; Davidson, Staplehurst & Dalzicl 2006; Staplehurst, Davidson & 
Dalzicl 2008). Moreover, a symmetry breaking between cyclonic and anticyclonic vortices 
is observed, both in forced (Hopfinger et al. 1982) and decaying (Bartello, Metais & 
Lesieur 1994; Smith & Waleffe 1999) turbulence. This symmetry breaking has received 
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considerable interest in recent years. It has been quantified in terms of the vorticity 
skewness, = (wf ) / (ui 2 ) 3 ^ 2 (where io z is the vorticity component along the rotation 
axis), which is found to be positive for Ro ~ 1 (Bartello et al. 1994; Morize et al. 2005; 
Rupper-Felsot et al. 2005; Bourouiba & Bartello 2007; Bokhoven et al. 2008; Staplehurst 
et al. 2008). In decaying rotating turbulence, starting from initial conditions such that 
Ro S> 1, a power-law growth has been observed in the form S u ~ (Sl£)°- 6±al by Morize 
et al. (2005). The fact that the time appears through the non-dimensional combination 
fit actually suggests a build-up of vorticity skewness acting on the linear time scale f2 _1 . 

Several explanations have been proposed for the cyclone-anticyclone asymmetry growth, 
although none provides a complete explanation of the experimental data. First, in a ro- 
tating frame, the vortex stretching term for the axial vorticity, (2Q + uj z )du z /dz, is larger 
for cyclonic than for anticyclonic vorticity. Gence & Frick (2001) have shown that, for 
isotropic turbulence suddenly subjected to a background rotation, S u grows linearly at 
short time, i.e. for t <C f2 -1 . At larger time the growth is expected to be slower, because 
the vertical strain du z /dz is reduced by the rotation. Second, anticyclonic vortices are 
more prone to centrifugal instabilities. This effect can be readily shown for idealized 
axisymmetric vortices, for which the generalized Rayleigh criterion in a rotating frame 
(Kloosterziel & van Heijst 1991), <fr(r) — 2(oj z + 2fl)(ug/r + f2), is more likely to become 
negative for anticyclonic vorticity. Sreenivasan & Davidson (2008) have shown, using 
a model of axisymmetric vortex patches, that cyclonic vortices first develop columnar 
structures, while anticyclonic vortices become centrifugally unstable. 

More surprisingly, a decrease of S u has been reported at large times (smaller Ro) by 
Morize et al. (2005, 2006b), and later confirmed by Bokhoven et al. (2008). There is no 
general agreement concerning this unexpected return to vorticity symmetry. Several two- 
dimensional mechanisms, such as cyclone merging or diffusion, may induce a reduction 
of this quantity. The role of the confinement was also suggested by Morize et al. (2005), 
nameling the non-linear Ekman pumping on the rigid walls. However it was questioned 
by Bokhoven et al. (2008), who also reported a decrease of 5 W at large times, but in a 
numerical simulation with periodic boundary conditions, hence with no Ekman pumping. 
The non-monotonic time evolution of the vorticity skewness is confirmed by the present 
experiment, although the present boundary conditions significantly differs from those of 
Morize et al. (2005) . The present results suggest another contribution for this decrease at 
large times: As time proceeds, the vertical velocity, initiated by the 3D initial conditions, 
forms vertically coherent layers, either ascending or descending, passively transported by 
the large-scale quasi-two-dimensional flow. The straining of those layers by the horizontal 
motion produces smaller scales, in a process similar to the enstrophy cascade in two- 
dimensional turbulence. This mechanism reinforces the vertical shear, making those layers 
prone to inertial instabilities, producing small-scale horizontal vorticity, which in turn 
produces either cyclonic or anticyclonic vertical vorticity, resulting in a reduction of S u at 
large times. This re-injection of symmetric vorticity fluctuations at small scale is thought 
to be a generic mechanism in decaying rotating turbulence, provided the initial state 
contains a significant amount of vertical velocity, which is the case for an initial 3D 
isotropic turbulence. 

The paper is organised as follows. Section [2] describes the experimental setup and the 
PIV measurements. In Sec. [3] the large scale flows and the procedure to extract the 
turbulent flow component are detailed. The influence of the background rotation on the 
energy decay and the time evolution of the non-dimensional numbers are presented in 
Sec. |4j The anisotropy growth and the formation of the vertical layers are characterised 
in Sec. [5] The structure and dynamics of the vertical vorticity field is described in Sec.[6j 
with emphasis given on the cyclone-anticyclone asymmetry growth and its decay induced 
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by the shear instability of the vertical layers. Finally, Sec. [7] summarizes the different 
regimes observed during the decay. 



2. Experimental set-up and procedure 

2.1. Experimental apparatus 

The experimental setup, sketched in figure [TJ consists in a 13 m x 4 m water channel 
mounted on the 'Coriolis' rotating platform. Details about the rotating platform may be 
found in Praud et al. (2005, 2006) and Praud & Fincham (2005), and only the features 
specific to the present experiments arc described here. The whole tank is filled with 
water to avoid leaks from the channel, with a depth of h = 1 m. One experiment without 
rotation, and 3 experiments with rotation periods of T = 30, 60 and 120 s, have been 
carried out (see table [TJ. The angular velocity = 2n/T is constant within a precision of 
Ail/ft < 10~ 4 . The parabolic deformation of the surface height induced by the rotation 
is 0.3 cm (resp. 4.5 cm) for the lowest (resp. highest) rotation rate. 

The turbulence is generated by horizontally translating a vertical grid at a constant 
velocity V g = 30 cm s _1 along the channel over a distance of L x = 9.1 m (sec Movie 
1). The grid is made of square bars of 30 mm, with a mesh size of M — 140 mm and 
a solidity ratio of 0.38. The grid cross section is 4 m x 1 m = 28M x 7M, ensuring 
negligible lateral confinement effects. It is hung from a carriage moving above the free 
surface, and its displacement is driven through a computer-controlled DC motor, ensuring 
good repeatability of the experiments. The velocity of the grid increases linearly from 
to V g and, remains constant in the central part, and decreases linearly back to zero 
at the end of the channel. The acceleration and deceleration are chosen to minimize the 
generation of gravity waves by accumulation of fluid between the grid and the endwalls 
of the channel. The time at which the grid reach the center of the channel, where the 
measurements are performed, defines the origin t = 0. The streamwise, spanwise and 
vertical axis are noted e x , e y and e z respectively. 

The initial condition of an experiment is characterized by two nondimensional param- 
eters, the Reynolds and Rossby numbers based on the grid velocity and grid mesh, 

Re n = — - — -, R°g = 



20M' 



where v is the water kinematic viscosity. The grid Reynolds number is constant for all 
the experiments, Re g — 4.20 x 10 4 , while the Rossby number lies in the range 5.1 — 20.4 
(table [JJ . Accordingly, the turbulent energy production in the near wake of the grid is 
expected to be weakly affected by the rotation, so the initial state can be considered as 
approximately isotropic. 

2.2. Particle Image Velocimetry measurements 

2.2.1. PIV setup 

A high-resolution particle image velocimetry system, based on a 14 bits 2048 x 2048 pix- 
els camera (PCO.2000), was used in these experiments. Water was seeded by Chemigum 
P83 particles, 250 /im in diameter, and illuminated by a 8W laser sheet of thickness 1 cm. 
The laser sheet was transmitted to the measurement area through an immersed mirror 
at 45° at one end of the channel. 

Two types of measurements were performed: 

(a) in a centered square area of 1.3 m x 1.3 m in the horizontal plane (e x ,e y ) at 
mid-height (z = 0.5 m). The camera is located 4 m above the horizontal laser sheet (CI 
in figure [JJ , and the area is imaged through the free surface (see Sec. 
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Figure 1. Side view (a) and top view (b) of the experimental setup. The grid is translated from 
left to right, along e x . The angular velocity is tt — Qe z , with fl > (anticlockwise rotation). 
The PIV camera is located either at CI or C2, for measurements in the horizontal and vertical 
planes, respectively. The dashed squares show the corresponding imaged areas. 



(b) in a 1.1 m x lm area in the vertical plane (e x ,e z ) in the middle of the channel. 
The plane is imaged through a window in the lateral wall (C2 in figure so that the 
measurements are not affected by free surface disturbances. 

Up to six decay experiments of 1 hour (7700 grid timescales M/V g ) have been carried 
out for each rotation rate and, for each decay, 400 image pairs are recorded. Since the 
characteristic velocity decreases in time, the delay between the two successive images of a 
pair is made to gradually increase during the acquisition sequence, from 125 ms to 2 s, so 
that the typical particles displacement remains approximately constant throughout the 
decay. The time delay between image pairs is also gradually increased during the decay, 
from 2 s to 20 s. The results are ensemble-averaged over the N r = 6 realizations in the 
horizontal plane, and N r — 4 in the vertical plane. Although this is enough to achieve 
statistical convergence at small times, when the correlation length is significantly smaller 
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Table 1. Flow parameters. The symbols are used in the following figures. 



than the imaged area, the convergence becomes questionable at large times, when the 
imaged area contains in average 1 large-scale structure or less. 

2.2.2. PIV resolution 

The PIV computations have been performed using the software Davi^ .and the statis- 
tical analysis of the velocity fields using the PivMat toolbox under MatlaqJ. Interrogation 
windows of size 32 x 32 pixels, with an overlap of 16 pixels, were used. The final velocity 
fields are defined on a 128 x 128 grid. The spatial resolution is Ax = 10 mm, which is of 
the order of the laser sheet thickness. With a particle displacement resolution of 0.1 pixel, 
achieved using a classical subpixel interpolation scheme for the correlation function, a 
velocity signal-to-noise ratio of 2 x 10~ 2 is obtained. 

Due to the moderate spatial resolution of 10 mm, the velocity field inside the Ekman 
boundary layer, of thickness 5g = (v/Q) 1 / 2 ~ 2.2 — 4.4 mm (see Table [TJ, cannot be 
resolved. Assuming isotropy in the bulk of the flow, which is valid only in the non-rotating 
case or at small time, the smallest turbulent scale can be estimated by the Kolmogorov 
scale rj — (i^/e) 1 / 4 , where the dissipation rate e can be computed from the energy decay, 
e ~ — (3/2)d(u' x ) 2 /dt (energy decays are detailed in Sec. |4j. The scale rj is of order of 
0.4 mm ~ Ax/30 at t ~ 20M/V g , for all rotation rates, so that the smallest scales 
are not resolved at the beginning of the decay. Accordingly, the measured vorticity at 
scale Ax underestimates the actual one. Vorticity measurements become reliable when 
rj > Ax/5, which is satisfied for t > 250 M /V g only. At the end of the decay, r\ is of order 
of 8 mm for Q — 0, in which case the vorticity can be accurately computed from the PIV 
measurements . 

2.3. Free surface disturbances 

When the flow is imaged from above, a significant uncertainty originates from the refrac- 
tion through disturbances at the free surface. These perturbations induce an additional 
apparent particle displacement, 5xfsd = (1 — l/n w ) (h/2) V/i, with n w the water refrac- 
tion index, V/i the surface gradient, and h/2 the path length of the refracted light rays 
(see e.g. Moisy, Rabaud & Salsac 2009). 

Two sources of free surface disturbances of small wavelengths have been identified: 
perturbations due to the wake of the grid, which prevent the measurements during typi- 
cally 10 s (~ 20M/V g ) after the grid translation, and surface waves continuously excited 
by residual mechanical vibrations from the platform structure, in the range 5-10 Hz 
(wavelength of 10-30 cm), which may be non negligible for the highest rotation rate. The 
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Figure 2. Horizontal and vertical snapshots of the velocity fields taken at t = 360 s ~ 770M/V g 
after the grid translation, without rotation (left figures) and with rotation at Q — 0.10 rad s _1 . 
The imaged area is 1 x 1 m. The color shows the corresponding normal component of the 
vorticity, ui z (x, y) and uj v (x, z). For the rotating cases, this time corresponds to 6 tank rotations. 
The color palette for the vorticity ranges from —0.1 to 0.1 rad s" 1 . Note that mean flow in 
the direction of the grid velocity (along e x , which is present for Q = but inhibited by the 
background rotation for 7^ 0. 

resulting velocity contamination is of order of |ufsd| = |<5xfsd|/<^, with St the inter- 
frame time. The disturbance amplitude was significantly reduced when a thin layer of 
impurities (dust) was present at the water surface. In practice, the velocity contamina- 
tion originating from those free surface disturbances can be neglected when computing 
quantities based on the velocity field itself, or on the antisymmetric part of the velocity 
derivative tensor (e.g. vorticity), but may significantly affect its symmetric part (e.g. 
horizontal divergence). 

3. Large scale properties of the flow 

3.1. Flow visualisations 

First insight into the influence of the background rotation on the turbulence decay may 
be obtained by comparing the horizontal and vertical vorticity fields shown in figure [H 
in the non-rotating case (left) and in an experiment rotating at O = 0.10 rad s _1 (right). 
Those snapshots are obtained 360 s after the grid translation (770M/V g ). At this time, 
the turbulent Reynolds number (defined in Sec. I4.4[) is 400 and 700 for the non-rotating 
and rotating cases respectively, and the Rossby number for the rotating case is 0.06. 

While the vorticity fields u z and u> v for the non-rotating cases are similar in the two 
measurement planes, as expected for approximately isotropic turbulence, they strongly 
differ in the rotating case. Movies 2 and 3 (see supplementary matarial) of lo z and u> y 
clearly show the two essential features of the turbulenc decay in the rotating frame, 
namely the anisotropy growth and the cyclone-anticyclone asymmetry. 

The vertical vorticity co z shows strong large-scale vortices, mostly cyclonic (in red), 
surrounded by shear layers. In the vertical plane, the spanwise vorticity uo y shows elon- 
gated structures of alternating sign, originating from layers of ascending and descending 
fluid. The dominant contribution of uj y comes from the vertical shear, du z /dx, except 
near the top and bottom boundary layers where the horizontal shear du x /dz is dominant. 
Sections [5] and [6] will be devoted to a more detailed characterisation of the structure of 
the flow. 
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Figure 3. Time evolution of the ensemble and spatially averaged streamwise ( — ) and spanwise 
(- -) velocity components, for Q — 0.1 rad s _1 (rotation period T = 60 s). The fast oscilla- 
tion, of period Taw — 7.3 s, is a longitudinal gravity wave, and the slow oscillation, of period 
Tiw = T/2 — 30 s, is an anticyclonic inertial wave. 



3.2. Large scale flows 

Translating a grid in a closed volume is ideally designed to produce homogeneous turbu- 
lence with zero mean flow. However, repeatable flow features are found over successive 
realisations, so a careful separation between mean and turbulent flows is necessary to 
analyse the present experiments. 

The main repeatable features of the flow generated by the grid translation can be 
inferred from figure [3l where the spatially- and ensemble-averaged streamwise, (U x ), and 
spanwise, (U y ), velocity components are shown as a function of time, for fl = 0.1 rad s _1 
(the averaging procedure is detailed in Sec. l3.3[) . A large scale circulation (LSC) is clearly 
visible, traced here as a systematic positive, slowly decaying, streamwise velocity. In 
addition, two oscillatory flows, namely a gravity wave (GW) and an inertial wave (IW), 
are found to superimpose to the LSC. 

Gravity wave (GW): As the grid is translated along the channel, it pushes a significant 
amount of water near the endwall, which initiates a fast longitudinal gravity wave. This 
sloshing mode is present both in the rotating and non-rotating cases. Its period Tq W is 
given by A/c, where the wavelength A is twice the channel length and the phase velocity 
in the shallow-water approximation is c = (gh) 1 ^ 2 ~ 3.13 m s _1 , yielding Tqw — 7.3 s 
(see figure [3]) ■ 

Large scale circulation (LSC): This mean flow is responsible for the positive mean 
streamwise velocity at mid-height in figure [3j of initial amplitude of order ~ 2 x 10~ 2 V g . 
It originates from the boundary condition asymmetry between the solid boundary at 
z = and the free surface at z — h, yielding a significant residual horizontal shear 
d(U x )/dz > 0, of initial amplitude ~ 0.03V g /h ~ 10~ 2 s -1 . This shear flow is indeed 
visible in figure |^c) for SI = 0. This positive streamwise flow is also present in the 
depth-averaged velocity measured in the vertical plane, and must be compensated by a 
horizontal recirculation flow along the lateral walls. 

Lnertial wave (IW): When rotation is present, the mean horizontal shear produced by 
the grid further excites an inertial wave, of period half the rotation period of the tank 
Tiw = T/2 (Thy = 30 s in figure[3]). Visualisations in the vertical plane show that the IW 
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can be approximately described as the oscillating shearing motion of horizontal layers, 
of thickness h/2, so that the associated wavenumber is vertical. In the horizontal plane, 
the signature of this IW is a uniform anticyclonic oscillation, visible by the phase shift 
of 7r/2 between the mean velocity components (U x ) and (U y ). Its amplitude is of the 
order of Uxw — lO^ 2 !^ (resp. 5 x 10~ 4 Vg) at the beginning (resp. end) of the decay. The 
gyration radius associated to this wave is f/rw^iw /^tt — 3 cm (resp. 1 mm). 

In the non-rotating case, the mean horizontal shear persists over large times. Although 
very weak, it acts as a source of turbulence, which is found to affect the energy decay 
at large times (Sec. 2]). On the other hand, when rotation is present, this mean shear 
oscillates at a frequency 2f2 which is about 100 times larger than d(U x )/dz, yielding a 
strong suppression of the turbulence production. A nice consequence of this effect for the 
present experiment is that translating a grid in a rotating frame produces a turbulence 
which is much closer to homogeneity than in a non-rotating frame. 

Finally, we note that the disturbance height associated to both the GW and IW flows, 
estimated as £ = {h/gY^U where U is the wave velocity, is less than 0.6 mm, yielding 
negligible surface slopes. As a consequence, contrarily to the free surface disturbances 
of small wavelength discussed in § 12. 3( no significant optical distortion is expected here 
from those large scale waves. 

3.3. Large scale flows subtraction 

To summarize, the measured velocity field U can be written as the sum of 5 contributions: 
the 3 large scale flows (LSC, GW and IW), the turbulent field of interest u, and the 
apparent velocity field induced by the free surface disturbances ufsd (present only at 
early time for the measurements in the horizontal plane) . The three large scale flows are 
simply described by their amplitude, U^sc, Uqw an d Urw, which are slowly decaying 
in time. Keeping only the dominant spatial dependences of those 5 contributions, the 
measured velocity can be approximated as 



TjW (x,y,z,t) ~ U LSC (z,t)e x + U GW (t) cos (jr^J 



+Utw(z,€) 



( 2-Kt \ ( 2vrf \ 

-u {n) { Xl y,z,t) + uf^ D (x,y,z,t), (3.1) 



where n is the realization number (the phase origin of the GW and IW flows are not 
considered for simplicity) . Within this approximation, the three large scale flows are uni- 
form translations in the horizontal plane. As a consequence, providing that the turbulent 
scale is significantly smaller than the imaged size, these large scale flows can be readily 
subtracted from the measured velocity fields. For the measurements in the horizontal 
plane, neglecting ufsd, one has: 

uW ( x , y , t) = TJW (x, y,t)- {Ujri (x, y,t)) x , y , B , (3.2) 

with a = x,y. Here the brackets (-} x .y,e denote both ensemble and spatial averages. For 
a field A^ n '(r,t) the ensemble average (-) e is defined as 



1 N r 



where n = l..N r is the realisation (N r — 4 and 6 for the measurements in the vertical 
and horizontal plane respectively). The spatial average along the direction x a is defined 
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Figure 4. Total and turbulent kinetic energy (streamwise variance), (a), f2 = 0. (b), 
= 0.10 rad s" 1 (T = 60 s). The oscillations in the rotating case correspond to IW flow, 
of period Trw = T/2 = 30 s. 
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where L a is the image size in the e a direction. Similarly, for the measurements in the 
vertical plane, the average {-) x ,z,e is computed. In order to reduce the statistical noise 
owing to the limited number of realisations, a temporal smoothing is also performed, on a 
range [ti,i2] such that t% — t\ < (t% +t2)/20. When there is no ambiguity, single brackets 
(•) denotes in the following both the ensemble and spatial averages and the temporal 
smoothing. Finally, the root mean square (rms) is noted A' = [A 2 ) 1 / 2 . 

It must be noted that, at the end of the decay, when the scale of motion becomes 
of the order or even larger than the imaged size, distinguishing the large scale and the 
turbulent contributions becomes ambiguous, and the subtraction may underestimate the 
actual turbulent energy. 

4. Energy and integral scales 

4.1. Energy decay 

The time evolution of the streamwise velocity variance for the total flow, (U 2 ), and the 
turbulent flow, (u 2 ) — ((U x — (U x )) 2 ), are shown in Figured] In the absence of rotation 
(figure 2k), the energy of the LSC flow clearly dominates the total energy, by a factor 
up to 10 for t ~ 1000M/Vg. On the other hand, when rotation is present (figure HJa) , 
the turbulent energy is very close to the total energy, confirming that the mean flow is 
significantly reduced in the presence of rotation. 

Once the mean flow is subtracted, the turbulent energy in the non-rotating case decays 
as t~ n up to t ~ i00M/Vg, with n ~ 1.22 ±0.05 (figures H and [SJa,) . This decay exponent 
is very close to the Saffman (1967) prediction n = 6/5 for unbounded turbulence. The 
streamwise variance, {u 2 ,), is approximately 1.4 times larger than the 2 spanwise variances 
(v,y) and (u 2 ), reflecting the usual residual anisotropy of grid turbulence (Comte-Bellot 
& Corrsin, 1965). For t > 400Af/V g , however, the shallower decay originates from the 
turbulence production by the mean vertical shear, which is a specific feature of the non- 
rotating case. This transition time, noted t s hoar in figure [5]( a) , is indeed of the order of 
the shear timescale, (dUx/dz) -1 ~ 250 s ~ 540M/Vg. The ordering of the 3 velocity 



Anisotropy and cyclone- anticyclone asymmetry in decaying rotating turbulence 11 



: 1 1 1 


<u*> 






r :: ' : * ;: x 


— <u*> 




~ vv 




1 1 






Figure 5. Time evolution of the variance of the 3 velocity components. The variances (ujj) 
and (Uy) are computed from the horizontal PIV fields (camera CI), and {u 2 z ) from the vertical 
PIV fields (camera C2), for non-simultaneous experiments, (a), Q, — 0. The vertical tick i a hear 
indicates the time after which the turbulent energy production by the residual mean shear 
becomes significant, (b), Q = 0.05 rad s . The vertical tick t* indicates the transition between 
the i -6 / 5 isotropic decay and the t~ 3 ^ 5 decay affected by the rotation. 



variances for t t s hcar, > (u 2 ) > {u 2 ), actually confirms the shear-dominated 

nature of the turbulence in the non-rotating case at large times (Tavoularis & Karnik, 
1989). 

In non-dimensional form, the decay law of the streamwise variance (neglecting possible 
time origin shift) writes 

-6/5 

n^) (4.i) 



V 2 
<j 



A best fit yields a decay coefficient A ~ 0.045 ± 0.005 for t < i s hcar, a value in good 
agreement with the literature for grid turbulence (see, e.g., the review by Mohamed and 
LaRue, 1990). This low value reflects the weak efficiency of a grid to produce turbulent 
fluctuations. This demonstrates that, in spite of the significant large scale flow generated 
by the present forcing, the decay of the turbulent kinetic energy is close to that of classical 
grid turbulence, suggesting a negligible coupling between the mean flow and the small 
scale turbulence, at least for t < t s hear- 

The time evolution of the 3 velocity variances in the rotating case are shown in fig- 
ure [5jb) for = 0.05 rad s _1 . At early time, the 3 curves are very close to the reference 
case O = (figure [5^) , confirming that the rotation has no measurable effect at large 
Rossby numbers. After a crossover time t* ~ 100M/V g , the decay of the 2 horizontal 
variances (u 2 ) and (u 2 ) become shallower, showing a clear reduction of the energy decay 
by the rotation. On the other hand, the vertical variance (u 2 ) first follows the horizontal 
variance short after the crossover time t* , but sharply decreases soon after, reflecting a 
growth of anisotropy. Since here the turbulence production by the mean shear is essen- 
tially suppressed by the background rotation, this departure from the t -6 / 5 decay and 
the resulting anisotropy growth can now be interpreted as a pure effect of the rotation. 
The anisotropy of the flow will be characterized in more details in Sec. and we focus 
here on the behaviour of the horizontal variances only. 



4.2. Crossover between the two decay regimes 
In order to characterize the influence of the rotation on the transition time t* 



the de- 



cays of the streamwise velocity variance (u 2 ) are compared in figure [6] for the 4 sets 
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Figure 6. Time evolution of the streamwise velocity variance for the non-rotating and 

the 3 rotating experiments. The solid line shows A(tV g /M)~ 6 ' 5 and the dashed lines show 
AnRog 3/5 (tV g /M)~ 3/5 . The transition between the non-rotating (t~ 6/5 ) and rotating (t~ 3/5 ) 
decay laws occurs at t* , indicated by the 3 vertical ticks for each rotation rate. For t > tshear the 
turbulent energy production by the residual mean shear becomes significant in the non-rotating 
case. 



of experiments. The crossover time t* decreases from 100 to approximately 30M/V g as 
Q is increased, corresponding to approximately 0.4 tank rotation. The small value of t* 
found for the highest rotation rate (fi = 0.20 rad s _1 ) indicates that the turbulent energy 
production in the wake of the grid may be indeed already affected by the background 
rotation in this specific case (the grid Rossby number is Ro g = 5.1 only). The energy 
decay for this high Q is indeed particular, showing unexpected large fluctuations of the 
streamwise velocity variance. 

In the limit of large rotation rate, the energy decay can be modeled by assuming that 
the energy transfer rate scales as the linear timescale f2 _1 . Based on this argument, 
Squires et al. (1994) proposed, using dimensional analysis, the following asymptotic de- 
cay law, 

l,.2\ /*!/• \ -3/5 

AnRo-^r-^-) , (4.2) 



with Aq a non-dimensional constant. Although the elapsed time is moderate here, the 
decay curves in figure [6] are actually compatible with this shallower decay (|4.2[) . Fitting 
the data for Q, = 0.05 and 0.10 rad s^ 1 yields Aq. ~ 0.020 ± 0.005 (the data at fl = 
0.20 rad s" 1 being excluded for the reason given before). Accordingly, the crossover 
time t* between the non-rotating and the rotating decay laws is obtained by equating 
Eqs. gIE]) and (jO) . 



M \A n J 



5/3 



Ro g ~ (5± l)Ro g , (4.3) 



yielding the values 100, 50 and 25 for the 3 rotation rates, which reproduces correctly 
the observed t* in figured! at least for ft — 0.05 and 0.10 rad s _1 (see the vertical ticks). 
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Expressing this crossover time (|4.3j) in terms of the rotation rate is consistent with a 
transition occurring at fixed fraction of tank rotation, 

Qt* 

~ (5±1)/4tt~0.4±0.1. 

2tt 

This result indicates that, after only half a tank rotation, the turbulence decay is domi- 
nated by the rotation. 

It is remarkable that the transition between the two regimes i~ 6 / 5 and i" 3 / 5 is suffi- 
ciently sharp, so that the analysis of Squires et al. (1994) can be recovered to a correct 
degree of accuracy. A similar transition in the form t~ 10 / 7 — > f~ 5 / 7 , with again a factor 
2 between the non-rotating and the rotating decay exponents, has been observed in the 
recent simulation of Bokhaven et al. (2008), the discrepancy with the present exponents 
being probably associated to different energy content at small wavenumber in the exper- 
iment and the simulation. The steeper decay laws reported by Morize & Moisy (2006a), 
with a transition from t~ 2 to t^ 1 as f2 is increased, is probably an effect of the significant 
lateral confinement present in that experiment. 

4.3. Integral scales 

The integral scales in the horizontal plane are useful to define the instantaneous turbulent 
Reynolds and Rossby numbers. They are defined as 

/>oo 

L aa A t )= C aa ^(r,t)dr, (4.4) 
Jo 

from the two-point correlation function of the a velocity component along the [3 direction, 

n , A (u a (x,t)u a (x + re fj ,t)} 

C aa> p(r,t) = t 2~\ • l 4 - 5 ) 

\ u al 

Computing L aa ^ from experimental data using directly Eq. (|4.4[) is not possible, because 
of the limited number of independent realizations for separations r approaching the image 
size. To limit this uncertainty, the integral in Eq. (|4.4[) is truncated at a suitably defined 
decorrelation length r* , 

L aa .p = \ C aa ^(r)dr, with C aa j(r*) = C tr . (4.6) 
Jo 

A threshold value of Ctr = 0.2 was found to give reliable results throughout the decay. 
At early time, when a reliable estimate of the full integral (|4.4[) can be achieved, the 
truncated estimate is found to underestimate the actual one by a factor of about 1.3 to 
1.5. Although this definition systematically under-estimates the actual integral scales, 
the trends observed from these quantities are expected to represent the evolution of the 
true length scales. 

It must be noted that the presence of large scale flows discussed in § 13.21 would imply 
an unphysical increase of the horizontal correlations C aa ^{r), and consequently of the 
integral scales, if computed directly from the total measured velocity U. Accordingly, 
the subtraction of the large scale flows is critical for reliable measurement of the integral 
scales. 

The time evolution of the longitudinal integral scale, averaged over the two horizontal 
directions x and y (noted here 1 and 2 by convention), 

Lf = — (£11,1 + £22,2), 
is plotted in figure O and shows little influence of the background rotation. A best fit 
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Figure 7. Time evolution of the longitudinal integral scale in the horizontal plane, 
Lf — (in.i + 1/22, 2)/2, for the 4 series of experiments. The upper and lower dashed lines 
show the scalings t 2 ^ 5 and t 1 / 6 for reference. 

yields a power law L/(t) ~ £0.35±o.05 ; for t < 1000M/V S , for all rotation rates. The 
scatter at larger time is probably a consequence of the inadequate subtraction of the 
mean flow, which may occur when the size of the largest vortices becomes comparable 
to the imaged area. 

Dimensional analysis actually predicts different growth laws for Lf in the non-rotating 
and rotating cases (Squires et al. 1994), 

with B and Bq, non-dimensional constants. Surprisingly, although the i~ 6 / 5 — ► i~ 3 / 5 
transition at t — t* is evident in the energy decay curves (figure [6|) , there is no evidence 
for the equivalent t 2 / 5 — > i 1 / 5 transition for Lf in figure [7) Within the experimental 
uncertainty, a single power law t 2 / 5 actually provides a reasonable description for the 
growth of Lf both in the non-rotating and in the rotating cases. 

4.4. Instantaneous Reynolds and Rossby numbers 

The instantaneous Reynolds number, and the macro and micro Rossby numbers (Jacquin 
et al. 1990), are finally defined as 

Re(t) = <±l, Ro{t) = ^ Ro u (t) = ^. (4.8) 

The time evolution of those numbers are plotted in figure [5£a-b). After a short period of 
sharp decay similar to the non-rotating case, the Reynolds number in the rotating cases 
show a very weak decay in the range t* < t < 3000 Mj V g . Typical Reynolds numbers 
range from 500 to 1300 as the rotation rate is increased. 

In spite of the absence of clear transition in the growth law for the integral scale, the 
time evolution of Re(t) compares correctly with the forms expected for the non-rotating 
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Figure 8. (a), Reynolds number Re(t) = u' x Lf/v. (b), Micro- and Macro- Rossby numbers, for 
the three experiments with background rotation. Upper curves (black symbols): Ro^ = u' z /2Q. 
Lower curves (open symbols): Ro = u' x /2QLf. The horizontal dotted lines show the thresholds, 
Rou, = 1.8 and Ro — 0.25, with the corresponding transition times t* indicated by the vertical 
ticks [see Eq. (|4.3[) ]. The dashed lines show the scalings t _1 and i -1 / 2 for reference. 



and rotating cases - see Eq. (|iTT]) . P~2")) and l[3~7|) . 
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The micro-Rossby number Ro u , shown in figure [Hb), takes values about 10 times 
larger than Ro throughout the decay (note that Ro u may be underestimated at small 
times because of the limited PIV resolution). This moderate ratio indicates that the 
range between the large scales dominated by the rotation and the small scales is indeed 
limited for the Reynolds number of the present experiments. 

The joint decay and growth laws for the velocity and integral scale actually lead 
to a remarkably simple decay law for the macro Rossby number Ro(t). Combining 
again Eq. (|4.1|) . (|4.2[) and (|4.7|) shows that, for t < t* , the non-linear time scale r„; = 
Lf(t)/u' x (t) is simply proportional to the elapsed time t, with no dependence on the 
initial grid time scale M/V g . Since the Rossby number is given by Ro = (2ilT n i)^ 1 , it 
turns out to be only a function of the number of tank rotations, 



Ro(t) (x Ro g ( ^ ) oc (2fit) 



M 



(t « t*), 



tv„ 



-1/2 



Ro{t) oc Rol' 2 cx(2r>i)- 1 / 2 (t»f). 

At the transition t = t* , which is reached after a fixed number of rotations, the Rossby 
number is indeed found approximately constant, Ro(t*) ~ 0.25 (see the vertical ticks at 
t* in figure [HJd) . This value is in correct agreement with the transitional Rossby numbers 
reported by Hopfinger et al. (1982) and Staplehurst et al. (2008). 

No such simple law applies for the micro-Rossby number Ro u (t) . For t <C t* , assuming 
again isotropic turbulence, the decay law of Ro^ can be inferred from the relation between 
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Figure 9. Sequence of 6 snapshots of the velocity and spanwise vorticity ui y in the vertical plane 
(x, z) for — 0.20 rad.s -1 . The imaged area is 1 m x 1 m. The grid is translated from left to 
right, and the time origin t — is defined as the grid goes through the center of the imaged 
area. The color range is normalized by the rms uj' y = (cuy)w^, e computed for each time. 



the vorticity rms, the velocity rms and the dissipation rate, 



Idu' 2 



3 du'; 



2 dt 2 dt 
Combining the isotropic decay law (|4. 1(1 with Eq. (|4.9p yields 



(4.9) 



-11/10 



(t « **)> 



with a scaling exponent very close to that of Ro(t). Evaluating Ro^ at the transition 
t ~ t*, using Eq. finally yields 



Ro^t* 



(5±1) 



- 11 l w Ay 2 Re l J^Ro- 1 ' w . 



Accordingly, no strictly constant micro-Rossby number is expected at the transition, al- 
though the dependence with the rotation rate, as fi 1 / 10 , is very weak (Q is varied by a 
factor of 4 only in the present experiment) . As shown by the vertical ticks in figure \S[b) , 
Rou takes values which actually turn out to be approximately constant at the transition, 
Roui{t*) ~ 1.8. Interestingly, this value is close to the empirical threshold reported by 
Morize et al. (2005), below which the energy spectrum and the velocity derivative skew- 
ness were found to depart from the classical Kolmogorov predictions (computation of 
the macro- Rossby number based on the integral scale were not possible in that study). 
Although the macro-Rossby number is probably a more relevant parameter to describe 
this transition, the similar Reynolds number of the two experiments explains the similar 
values of Ro^, found at the transition. 



5. Dynamics of the anisotropy 

5.1. Visualisation of the vertical layers 

We now focus on the growth of anisotropy in the vertical plane (x, z). Figure [9] shows a 
sequence of 6 snapshots of the velocity field and spanwise vorticity uj v after the transition 
t > t*, for £1 = 0.20 rad s _1 (see also the Movie 3). The anisotropy can be visually 
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Figure 10. (a), Time evolution of the vertical velocity variance (iff), for the non-rotating and 
the 3 rotating experiments, (b), Isotropy factor u' z /u' x for the 3 rotating experiments, as a 
function of number of tank rotation. The vertical line indicates the transition between the £ -6 / 5 
and the t -3 ^ 5 decay regimes at Qt* /2n ~ 0.4. 



detected from the first snapshot, and the presence of vertical layers of ascending or 
descending fluid becomes evident after 8 tank rotations (figure [9}:). Although those layers 
are difficult to infer from the velocity field itself, because of the superimposed strong 
horizontal flow, they clearly appear through the surrounding layers of nearly constant tu y 
of alternate sign. Those layers of vertical velocity are consistent with a trend towards a 
three-component two-dimensional flow, with zero vertical variations of the velocity field, 
du/dz = 0, but non-zero vertical velocity u z originating from the initial conditions. 

As time proceeds, the vertical layers become thinner and more vertically coherent 
(note that since only the intersection of the layers with the measurement plane can 
be visualized, the apparent thickness may overestimate the actual one). At large time 
(figure E]3,f), although those layers are nearly coherent from the bottom wall up to the 
free surface, they are not strictly vertical, but rather show wavy disturbances. Those 
disturbances have amplitude and characteristic vertical size of the order of the layer 
thickness, suggesting the occurrence of a shear instability. We will examine in Sec. [6] the 
consequence of this instability on the dynamics and statistics of the vertical vorticity 
field. 

5.2. Decay of the vertical velocity and anisotropy growth 

The time evolution of the vertical velocity variance v! z = (Ug) x ,z,e and the isotropy ratio 
u' z /u' x , plotted in figure llOT a-b). show a complex and non-trivial behaviour. Here, the 
spatial average is computed only in the core of the flow, excluding layers of thickness 
O.l/i near the bottom wall and the free surface. 

Similarly to the horizontal variance (see figure [H]), the vertical variance for the rotating 
cases first departs from the reference curve i" 6 / 5 of the non-rotating case, and follows a 
shallower decay which is compatible again with a i~ 3 / 5 law, at least in an intermediate 
temporal range. Although the i -6 / 5 — > t~ 3 / 5 transition is not as sharp as for the horizon- 
tal variance, perhaps because of the limited statistics achieved for the measurements in 
the vertical plane, the transition time t* is compatible with the one determined for (u x ), 
corresponding to fit* /2tt ~ 0.4 tank rotation. Since both horizontal and vertical com- 
ponents follow the same decay law t~ 3 / 5 short after t*, the flow remains approximately 
isotropic, as shown by the constant ratio u' z /u' x ~ 1 ± 0.1 in figure llOT b) ■ However, this 
isotropy holds only during a restricted time range, until Qt/2ir ~ 2 tank rotations. After 
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Figure 11. Time evolution of the Ekman ratio pE = 2u' z / 8eoj' z . 



this time, corresponding to approximately t ~ (400 — 800)M/V^ in figure [TW a). the ver- 
tical variance follows a significantly sharper decay, whereas the horizontal variance still 
decays as i -3 / 5 , yielding a growing anisotropy consistent with the flow visualisations in 
the vertical plane. 

Although the formation of vertical structures is evident in the spanwise vorticity field 
uj y (figure [9]) , the anisotropy remains however moderate when expressed in terms of the 
velocity variances, as shown in figure fTOT b). The ratio u' z /u' x reaches a weak minimum 
between 0.6 and 0.4 only, after 10 — 30 tank rotations, depending on the rotation rate. 
Interestingly, the ratio u' z /u' x for the different rotation rates collapse in the anisotropy 
growth regime when plotted as a function of the number of tank rotations Qt/27r, sug- 
gesting here again an essentially linear mechanism for the building of this anisotropy. 

A remarkable feature of figure fTDTb') is the reverse trend u' z /u' x — ► 1 observed at large 
time. A similar behaviour is obtained for the Reynolds stress anisotropy in the numerical 
simulations of Morinishi, Nakabayashi & Ren (2001). This apparent return to isotropy, 
which breaks the scaling with the linear timescale J! -1 , is associated to the flattening of 
the decay of (u z ) at large time, visible in figure [TUT a) . However, the corresponding flow 
fields at large time (figures [9^,f ) are clearly not isotropic, but rather show wavy thin 
vertical layers. 

It must be noted that the magnitude of vertical velocity measured throughout the 
decays, even at large time, remains comfortably larger than the one expected for the 
Ekman pumping induced by the horizontal flow. According to the linear Ekman pumping 
theory, a quasi-2D field of vertical vorticity rms w' z would lead to a characteristic vertical 
velocity rms of u' z = 5ew' z /2 f° r z — &e (Greenspan 1968), where Se is the Ekman 
layer thickness (see Sec. I2.2.2")) . The discrepancy between the actual u' z and the Ekman 
pumping estimate u'f may be therefore measured by the ratio 

Pe = i-^t- (5-1) 

Oe^ z 

Note that this ratio corresponds, within a numerical prefactor, to the ratio between the 
Taylor scale A = yl5u' z /cj' z and Se- Figure [TT] shows the time evolution of pe, computed 
from non-simultaneous measurements of uj' z and u' z in the horizontal and vertical mea- 
surements respectively. Values around 10 at early time, slightly increasing up to about 
20 at larger time, are obtained, confirming that effects of the Ekman pumping can be 
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Figure 12. Time evolution of the normalized integral scales L aoit p/M computed in the vertical 
plane, for fl = 0.10 rad s~ . By convention, Z/11,3 is taken equal the normalized tank height, 
h/M ~ 7 (upper dotted line), when the vertical correlation of the horizontal velocity, Cn,3, does 
not decrease below 0.2. 

neglected in the present experiments. Accordingly, the vertical fluctuations found at late 
time originate from the initial vertical fluctuations induced by the grid. 

It is remarkable that the time at which the ratio pe becomes constant corresponds 
approximately to the time t* at which the energy decay curves depart from the non- 
rotating reference curve. The initial increase of pe may be simply described from the 
decay laws of Sec. H From Eqs. dHJ) and flU]), it follows that p E (t) oc (2m) 1 / 2 for 
t < t* , in qualitative agreement with the data of figure I 111 at early time. At the transition 
time Vlt* /2ir ~ 0.4, the value of pe is therefore independent of Ro g . The fact that u' z 
and lo' z remains approximately proportional for t > t* is an indication that the dynamics 
of the vertical velocity follows, in average, the dynamics of the vorticity. This supports 
the picture that the vertical velocity behaves as a scalar field, passively adverted by the 
horizontal flow at large time. 

5.3. Integral scales in the vertical plane 

In order to relate the evolution of the vertical velocity variance to the formation, thinning 
and instability of the vertical layers, we now focus on the statistical geometry of those 
layers. For this, we have computed the 3 integral scales £11,3, £33,1 and £33,3 from 
the velocity fields in the vertical plane, using definitions (|4.5|) - (|4.6|) with a,f3 = 1,3. 
£11.3 characterizes the trends towards two-dimensionality, £33.3 the vertical coherence 
of the layers, and £33,1 the thickness of the layers. No reliable measurement of L\\ 1 
could be obtained from the vertical fields, because of the ambiguity of the subtraction 
of the horizontal LSC flow at large time: Large scale vortices having their axis out of 
the measurement plane produce strong horizontal velocity which, if subtracted, yield an 
unphysical decrease of £11,1. Here again, in Eq. (|4.5|) . the depth-average excludes lower 
and upper layers over a thickness of h/10, in order to avoid boundary effects. In the 
extreme case of an unbounded z-invariant 2D flow, the vertical correlations would be 
Caa,3{r) — 1, yielding £ QQ ,3 = 00. By convention, if C(r) does not decrease below the 
threshold Ct r — 0.2 (see Eq. (I4.6|l ). L aa _p is taken equal to the channel depth h. 

The typical time evolution of the 3 integral scales is shown in figure [l2l in the case 
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£1 = 0.10 rad s _1 . At short time, the longitudinal integral scale £33,3 is, as expected, larger 
than the two transverse ones (one has £33,3 = 2Ln,3 = 2L33.1 for isotropic turbulence). 
The most spectacular effect is the rapid growth of £11,3, characterising the vertical corre- 
lation of the horizontal velocity, which is a clear signature of the two-dimensionalisation 
of the large scales of the flow. This rapid growth leads to a saturation of the integral scale 
for tVg/M > 400. Although the data are too noisy to check the scaling of this saturation 
time with the rotation rate, it occurs roughly at a constant number of tank rotations as 
fi is varied, which again supports the idea of a linear mechanism: At early time, energy is 
contained at scale Ln.i ~ M, and at each tank rotation eddies grow vertically by wave 
propagation, so £11,3 increases by an amount of L\\ \ (Davidson et al. 2006). Accord- 
ingly, L11.3 is expected to reach the channel depth h after a number of tank rotations of 
order of h/M ~ 7, which is in qualitative agreement with the present observations. 

After the saturation of in, 3, the vertical correlation of the vertical velocity remains 
constant until the end of the experiment, with £33,3 ~ 2M, indicating a significant, al- 
though finite, vertical coherence of the ascending and descending plane jets. Note however 
that even strictly coherent thin layers of constant velocity would lead to finite integral 
scale £33,3, because the tilting of the layers by the oscillating shear of the IW flow strongly 
reduces the vertical correlation as the layers become thinner. 

A remarkable feature of figure [T2] is the sharp decrease of the horizontal correlation of 
the vertical velocity, described by £33,1, for tV g /M > 2000 (corresponding to f2i/27r ~ 
16 — 24 tank rotations), and its subsequent saturation to the very low value £33,1 ~ 
0.2Af ~ 30 mm at large time. As far as the integral scales are concerned, the turbulence 
becomes highly anisotropic in the final stage of the decay, showing a non-trivial ordering 
(see figures 171 and 172]) : 

£33,1 *C £33,3 ~ la 1,1 -C in, 3- (5-2) 
The low asymptotic value of £33,1 suggests that, in the final regime, the vertical velocity 
fluctuations have a well defined characteristic scale in the horizontal direction, i.e. there 
is no global vertical motion at scales larger than the thickness of the layers. This final 
value of £33,1, which provides an estimate for the average thickness of the layers, is found 
to slightly decrease, from 0.23M to 0.17M, as fi is increased, suggesting that the thinning 
of the layers induced by the horizontal straining motion due to the large scale vortices is 
stronger at high rotation rate. 

5.4. Thinning and instabilities of the vertical layers 

In spite of the strong anisotropy of the integral scales in the final stage (|5.2p . the velocity 
variances showed an unexpected return to isotropy, u' z /u' x — ► 1 ffigurefTUb) . This apparent 
return to isotropy is also present at small scales, as shown by the following two velocity 
gradient isotropy factors, 

^JL an d ylPy. (5.3) 

where j' z = {{du z / dz) 2 ) 1 / 2 is the rms of the vertical strain rate. The vertical strain rate 
plays an important role, as it is responsible for the stretching of the absolute vertical 
vorticity. In isotropic turbulence, both quantities are equal to 1 (the second equality 
follows from the classical isotropic relation e = 15^7^ 2 = vu: 12 = ivui' 2 , where e is the 
dissipation rate). For a two-dimensional flow with arbitrary vertical velocity, one has 
7 Z = 0, whereas oj y = is true only for a two-dimensional two-component flow. As a con- 
sequence, the two isotropy factors may be considered as signatures of the dimensionality 
and componentality of the small scales respectively. 

Figure [13] shows that the two velocity gradient isotropy factors first slowly decrease 
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Figure 13. Time evolution of the velocity gradient isotropy factors: (a), uj' y /uj' z ; (b), \/^j' z /u)' z . 
In each figure, the horizontal dashed line indicates the isotropic values, uj' y /uj' z — y/E-y' z /uj' z — 1. 




Figure 14. Sketch showing thinning and instability of the vertical layers of vertical velocity 
advected by the horizontal flow, (a), Nearly vertical layer of ascending fluid, u z > 0, strained 
in the vicinity of a large cyclone uj z > 0. (b), the layer becomes unstable, producing horizontal 
vortices, (c) These horizontal vortices produce random horizontal motion, and hence vertical 
vorticity of arbitrary sign. 

according to the linear timescale fi^ 1 , reaching a moderate minimum of about 0.5 at 
the largest rotation rate. The time of maximum anisotropy for these quantities is close 
to that for u' z /u' x , although slightly lower, and here again the scaling with the linear 
timescale is broken during the increase at large time. This plot shows that, in the final 
stage, the small scales are both three-dimensional and three-component, although not 
necessarily isotropic. 

Assuming that the vertical velocity u z behaves as a scalar field passively advected 
by the large scale horizontal flow provides a qualitative explanation for the increase of 
ijj' y /u)' z at large time. As sketched in figure fT4Ta). a layer of ascending fluid u z > in a 
horizontal strain field, for instance in the vicinity of a large vortex, is elongated along 
one direction and compressed along the other one, so it becomes thinner. In this process, 
u z is approximately conserved, but its horizontal gradient Vh u z increases, producing 
horizontal vorticity lo x and oj y which may reach, and even exceed, the vertical vorticity 
u> z . 

The increase of 7^, on the other hand, may be a consequence of the instability of those 
vertical plane jets. If the inertial timescale of the jets, (VhU 2 ) _1 ~ L 33i i/u' z , remains 
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Figure 15. Sequence of 6 snapshots of the velocity and vertical vorticity fields ui z measured in 
a horizontal plane (x,y) at mid-height for fl = 0.20 rad.s -1 . The imaged area is 1.3 m x 1.3 m, 
representing 4.6% of the tank section. The tank rotation is anticlockwise. Positive and negative 
vorticity indicate cyclones (in red) and anticyclones (in blue) respectively. The color range is 
normalized by the rms uj' z — (tJz)w!y,e computed for each time. 

smaller than the dissipation timescale, 1 /v, the jets may undergo shear instabilities, 
producing horizontal vortices (sketched in figure [l"4b ). as suggested by the visualisations 
in figure [9je-f ) . This condition is actually satisfied: The Reynolds number Rej based on 
those jets, defined as the ratio of the two timescales, writes 

R L^K = W^ e , 
v L f u' x 

where Re is the instantaneous Reynolds number defined in Eq. (|4.8p . With L33 i/L/ ~ 
0.1, u' z /u' x ~ 0.5, and Re ranging between 300 and 1300 in the final period of the decay 
(see figured]), one has Rej ~ 10— 10 2 , which is actually sufficient for a shear instability to 
develop. Little influence of the background rotation is expected on this shear instability, 
since the vertical velocity is unaffected by the Coriolis force (the resulting instability 
pattern, involving horizontal velocity, may however be affected by the rotation). The 
resulting wavy jets break the vertical invariance of u z , thus producing vertical strain j' z 
of the order of the vorticity u>' y , in agreement with figure fT3l 

All these results suggest that the flow structure at large time is fully three-dimensional 
and three-component, with isotropy factors (|5.3p close to that of 3D isotropic turbulence, 
although the large scales are highly anisotropic, as described by the ordering of the 
integral scales (|5.2[) . 

6. Cyclone- Anticyclone asymmetry 

6.1. Dynamics of the cyclones and anticyclones 

We finally turn to the structure of the vertical vorticity field in the rotating case, focusing 
on the issue of the cyclone-anticyclone asymmetry. The dynamics of the horizontal flow 
is illustrated by the 6 snapshots in figure [TBI for t > t* (see also the Movie 2). At the 
beginning of the decay, the vorticity field consists in small scale disordered fluctuations 
(Fig. [T5a.b). which gradually evolve into a complex set of tangled vortex sheets and 
vortices (Fig. 115b ). A set of well defined, nearly circular, cyclones gradually emerges and 
separates from the turbulent background (Fig. fT5bl). Anticyclones are also encountered, 
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Figure 16. Vorticity skewness Su, and (b) Vorticity flatness F u , as a function of the number 
of tank rotation Q.t/2-K. In (a), the dashed line shows the fit 0.45(f2i/27r) ' 7 . In (b), the dashed 
line indicates the value F& — 3 corresponding to a Gaussian field. 



but they are weaker and less compact than the cyclones. Those anticyclones do not appear 
to be specifically unstable compared to the cyclones, so that the asymmetry seems to 
originate essentially from the enhanced vortex stretching of cyclonic vorticity. 

At large time, the size of the cyclone grows, and merging of cyclones are frequently 
encountered, as illustrated by figure [ToT e) . At the same time, a background of small 
scale symmetrical vorticity fluctuations appears (Fig. [TBe.f). At the end of the decay, 
the flow essentially consists in those small-scale symmetric fluctuations, advected by the 
large-scale, mostly cyclonic, motions. 



6.2. Growth of vorticity skewness 

The gradual structuration of the vorticity field can be described by the vorticity skewness 
and flatness factors, 

(where the brackets denote horizontal and ensemble average), which are plotted in Fig. 1161 
Both S u and show a growth followed by a decrease. In the growth regime at early 
time, both quantities collapse when plotted as a function of the number of tank rotations 
Qt/2n, indicating that the build-up of vortex structures is a linear mechanism. Note that 
the residual oscillations visible at small times are associated to the large scale IW flow 
(see Sec. I3.2j) . of period flt/2w = 1/2. On the other hand, as for the isotropy factors, this 
rescaling with f2 _1 no longer holds during the decrease of and F w at large time. 

For t < t* , the vorticity skewness is essentially zero, within an uncertainty of ±10 -1 . 
For t > t* , it grows according to the power law 

nt^ 0J 



S u ~ 0.45 — 

which is in remarkable agreement with the one reported by Morize et al. (2005), both 
concerning the exponent and the numerical pre-factor. Although in both experiments 
turbulence is generated by the translation of a grid, the details of the geometry differ 
in a number of respects: here the grid velocity is normal to the rotation axis, and the 
aspect ratio is significantly lower (h/L y — 0.25 instead of 1.3). The collapse of for 
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the two experiments is a clear indication of a generic behaviour of this quantity (Morize 
et al. 2006b). 

The peak values of S u , between 1.5 and 3 for increasing 17, are significantly larger than 
those obtained in the experiments of Morize et al. (2005) and Staplehurst et al. (2008), 
and in the DNS of Bokhoven et al. (2008). The corresponding peaks of F^, between 8 
and 18, are much larger than usually measured in non-rotating turbulence at similar 
Reynolds number (see, e.g., Sreenivasan & Antonia, 1997), an indication of the strong 
concentration of vorticity in the core of the cyclones. 



6.3. The decay of vorticity skewness at large time 

The peaks of S u and F u occur at the time corresponding to the return to isotropy of the 
small scales, as identified in Sec. 15.41 For larger times, when the flow consists mostly in 
isolated large-scale cyclones, S u decreases back to 0, while F u recovers values around 4, 
similar to the begining of the decay. The important scatter in the decay is due to the 
very limited sampling: At large times, the number of strong vortices per unit of imaged 
area becomes less than 1, so the statistics become very sensitive to events of vortices 
entering or leaving the field. 

There is no general agreement concerning the decrease of S w at large time. It was 
attributed to confinement effects by Morize et al. (2005), namely the non- linear Ekman 
pumping of the cyclonic vortices (Zavala Sanson & van Heijst 2000). Linear pumping, 
which is valid in the limit of u> <§C 20, should imply identical decay of vorticity moments 
of both cyclonic and anticyclonic vorticity, and should have therefore no effect of So,. 
This suggestion was motivated by the fact that the time t max of maximum Suj was 
approximately following the Ekman time scale, t ma x — O.lh^Q)^ 1 / 2 . Fitting the times 
of maximum skewness for the present data would actually give similar values, although 
the spread of the maximum of S u prevents from a real check of the fi" 1 / 2 scaling. 
However, the fact that the Ekman pumping is shown to have no significant effect in the 
present experiment (see Sec. [SJ rules out this interpretation. The role of the confinement 
in the decay of S u is also questioned by the numerical data of Bokhoven et al. (2008), 
who have reported a decrease of S u at large times in a homogeneous turbulence with 
periodic boundary conditions, and hence with no Ekman pumping. Note that no decrease 
of S u was reported in numerical simulation of Bourouiba & Bartello (2007) and in the 
experiment of Staplehurst et al. (2008), probably because of their limited temporal range 
(the measurements of the latter were restricted to 3 tank rotations, whereas the decrease 
of Suj starts typically after 10 rotations here and in Morize et al. 2005). 

The explanation we propose here for the decrease of relies on the fact that vor- 
ticity is a small-scale quantity, whereas the cyclone-anticyclone asymmetry is defined 
by structures of increasing size as time proceeds. In particular, it is observed that the 
flow outside the cyclones is not smooth, but is made of small scale vorticity fluctuations, 
approximately symmetric. These vorticity fluctuations originate from the instabilities of 
the vertical shear layers strained by the horizontal large-scale flow described in Sec. 15.41 
The horizontal vortices resulting from this instability induce an horizontal straining flow 
(fig. 114b ). which may itself be unstable and produce vertical vorticity of random sign 
at small scale (fig. IT4"b ~). Accordingly, there is a possibility for the vorticity skewness to 
return to zero at large time, although the large scale field remains dominated by a set of 
big cyclones. 

The characteristic size of the vortical structures may be estimated from the horizontal 
integral scale of the vertical vorticity, L u = {L^ 31 + Lg 3i2 )/2. Here the integral scales for 
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Figure 17. Time evolution of the horizontal integral scale of vertical vorticity, L, 



the vorticity are defined similarly as those for the velocity, by modifying Eq. (|4. 6|) as 

L^=f ^M^ + ^t)) ^ (61) 
Jo K/ 

where the truncation scale r* has the same meaning as in Eq. (|4.6p . Figure [T7] shows the 
characteristic increase and decrease of L u when rotation is present, whereas it mono- 
tonically increases in the absence of rotation. The decrease occurs at tV g /M ~ 2000 for 
all rotation rates, which coincides with the sharp decrease of £33,1 ( figure fT2|) and the 
isotropy factors (figure [T"3f . and reaches values of order 0.2M ~ 30 mm, similar to those 
found for £33,1. This is a confirmation that the vertical vorticity field is dominated, at 
large time, by the small scale fluctuations induced by the instabilities of the vertical shear 
layers. 

The role of the symmetric small scale vorticity fluctuation in the decrease of S u must be 
addressed very carefully, because the vorticity field deduced from the PIV may obviously 
be affected by measurement noise. The curl of this PIV noise being essentially symmetric, 
it would yield a trivial reduction of S u . This is a delicate issue, because the scale L u 
of those vorticity fluctuation is only slightly larger than the PIV resolution (Sec. 12.2^2]) . 
However, the temporal coherence of those small scale fluctuations advected by the large 
scales is evident at the end of the movie 2, whereas PIV noise would generate vorticity 
patterns essentially uncorrelated in time. Those vortical fluctuations cannot originate 
from the apparent particle displacements associated to the free surface disturbances (see 
§. I2.3|) . which would yield a purely irrotational displacement field. The temporal coherence 
of the flucutations may also be inferred from the 4 snapshots in figure 1181 where sets of 
arbitrary chosen vorticity patterns (marked in dashed ellipses) can be easily tracked in 
time, confirming that they are essentially advected by the large-scale motions. 



6.4. Skewness of the filtered vorticity field 

In order to check more precisely the influence of the measurement noise on S u , we have 
computed S u from the filtered velocity u obtained by convolution of u with a Gaussian 
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Figure 18. Sequence of four w 2 -snapshots in the horizontal plane (x,y) at large time, show- 
ing the advection of the small-scale symmetric vorticity by the large scale horizontal motion 
(Q — 0.10 rad s _1 , t ~ 64M/V g ). Each image is separated by 20 s = T/3, and the field of view 
is 1.3 mx 1.3 m. The two ellipses track some arbitrary vorticity pattern in time. The angular 
velocity of the cyclonic structure in the bottom-left corner (white circle) is fl c ~ 0.009 rad s , 
corresponding to a local Rossby number of fi c /fl = 0.09. 




Figure 19. (a), Time evolution of the skewness of the filtered vorticity field, Sa, for different 
filter size r/, for Q, = 0.20 rad s . (b), Vorticity skewness as a function of the filter size, at 
times t — 8T (before the peak of S w ) and t = 90T (after the peak). 



kernel of size rf, 

u(x,y,t;rf)=JJ u(x', y',t) -j=- e ^ dx'dy'. 

In practice the integral is restricted to a square area of size 6rf. 

The time evolution of the skewness of the filtered vorticity, S^, is shown in figure [TW a) 
for various filter sizes r/, in the case ft — 0.20 rad s _1 . In the growth regime, increasing 
the filter leads to a decrease of S^, showing that the vorticity asymmetry is essentially 
contained at the smallest scale. In this situation, although the measured may un- 
derestimate the actual one because of the finite resolution of the PIV measurement, the 
vorticity skewness truly reflects the cyclone-anticyclone asymmetry at the smallest scales. 
On the other hand, after the peak of S u , the ordering of the curves is reversed, so that 
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filtering the vorticity field now increases the skewness, showing that now the asymmetry 
is carried by vortices at larger scales. This is consistent with figure [TH where a large 
cyclone containing small-scale vorticity fluctuations is shown (white circle). However, al- 
though the peak of S u is shifted to larger times, a decrease is still observed. The effect 
of the filtering is further illustrated in figure fT9l b). where S& monotonically decreases 
as rf is increased at t = 8T, whereas it shows a non-monotonic behavior at t = 90T. 
Interestingly, in this latter case, the maximum of Sq provides a rough estimate of the 
size of the vortices responsible for the cyclone-anticyclone asymmetry at this time. 

One may conclude that, although S u provides a suitable description of the vortex 
asymmetry in the growth regime, when the characteristic size of the vortices corresponds 
to the diffusive scale (the 'Kolmogorov scale' modified by the rotation), it is no longer 
appropriate as the vortex size grows at larger time, in which case S u is strongly reduced 
by the small scale symmetric vorticity fluctuations. This does not imply, however, that 
filtering at even larger scales would totally inhibit the decrease of S u , since other physical 
mechanisms, such as vortex pairing or diffusion, may also produce a reduction of S^,. 

7. Conclusion 

The present experiment aims to focus on the transition at Rossby number Ro ~ 
O(l) which occurs in the course of the decay of grid turbulence, initially approximately 
homogeneous and isotropic, in a rotating frame. Emphasis is given on the energy decay, 
anisotropy growth and asymmetry between cyclonic and anticyclonic vertical vorticity. 

The different steps of the decay can be summarised as follows: 

(a) During the first 0.4 tank rotation (between 25 and 100 M/V g ), the instantaneous 
Rossby number Ro is larger than 0.25 and turbulence is not affected by the background 
rotation. Once the large scale mean flow and waves are properly subtracted, the turbulent 
energy follows the classical decay law £~ 6 / 5 of isotropic unbounded turbulence. 

(b) After 0.4 tank rotation, Ro < 0.25 and the first effects of the rotation are triggered. 
Provided the grid Rossby number is large enough, the energy decay in this regime is 
found to be compatible with the f2 3 / 5 i~ 3 / 5 law proposed by Squires et al. (1994), which 
is based on the assumption of energy transfers governed by the linear time scale . The 
horizontal flow becomes strongly correlated along the vertical direction, with a saturation 
of the integral scale Ln,3 by the vertical confinement. Both the large-scale isotropy factor 
u' z /u' x and the small-scale ones uj' y /uj' z and ^/uj^ depart from their isotropic value, on 
the linear timescale although reaching only moderate value of about 0.5. A cyclone- 
anticyclone asymmetry develops by preferential vortex stretching of the cyclonic vorticity, 
and is well described by a power law growth of the vorticity skewness as S u oc (ftt) ' 7 , 
consistent with the previous findings of Morize et al. (2005). 

(c) Finally, for t > 2000Af/V g (corresponding to 10-30 tank rotations), the vertical 
layers of vertical velocity, essentially advected by the large-scale horizontal flow, becomes 
thinner and prone to shear instabilities. These instabilities produce small scale spanwise 
vorticity, resulting in an apparent return to isotropy, which breaks the scaling with the 
linear timescale. However, the flow structure remains strongly anisotropic, as revealed by 
the characteristic ordering of the integral scales. A remarkable consequence of the insta- 
bility of those vertical shear layers is that it re-injects horizontal velocity disturbances, 
and hence vertical vorticity with random sign, at small scales. Those symmetric vorticity 
fluctuations appear as a small-scale noise, which contributes significantly to the reduction 
of the vorticity skewness S u , although the large scale vortices still remain preferentially 
cyclonic. 

The mechanism presented here provides an explanation for the decrease of S u at large 
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time, which is observed here and in other configurations (Morize et al. 2005; Bokhoven 
et al. 2008). It must be noted that, although S u provides a suitable description of the 
vortex asymmetry in the regime (b), when the vortices responsible for the asymmetry are 
at small scale, it is no longer appropriate as the vortex size becomes larger, in the regime 
(c), when the flow consists in small-scale symmetric vorticity fluctuations superimposed 
to large scale quasi-horizontal motions. A more suitable statistical quantity, based for 
instance on the low-pass filtered vorticity or the transverse velocity increments, should 
provide a better description of the cyclone-anticyclone asymmetry in this regime. 

Finally, the present results suggest that the initial conditions have a critical importance 
in the asymptotic state of decaying rotating turbulence, and in particular in the cyclone- 
anticyclone asymmetry. For initial isotropic turbulence, as is approximately produced in 
the wake of a grid, the initial vertical fluctuations, which represents 1/3 of the initial 
turbulent kinetic energy, is temporarily stored by the horizontal quasi-2D motions and 
plays no role in its dynamics, except at large time when it is released (and dissipated) 
directly at small scale via the shear instabilities of the vertical layers, resulting in a short- 
cut of the energy cascade. The effects of this energy re-injection at small scale is probably 
enhanced in the experiments of Morize et al. (2005), in which the grid was translated 
vertically instead of horizontally, resulting in a stronger decrease of the vorticity skew- 
ness at large time (Morize et al. 2006b). On the other hand, a rather different situation is 
expected for decaying rotating turbulence starting from strictly two-dimensional initial 
conditions, as in the experiments of Longhetto et al. (2002) and Praud et al. (2006), in 
which no vertical velocity is produced by the translation of a rake instead of a grid. In this 
latter situation, the vorticity skewness should reach significantly larger values, although 
a decay at large time may still be observed but for other reasons, such as vortex pairing 
or diffusion. 
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